Influence of extended dynamics on phase transitions in a driven lattice gas 
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Monte Carlo simulations and dynamical mean-field approximations are performed to study the phase 
transition in a driven lattice gas with nearest-neighbor exclusion on a square lattice. A slight extension of 
the microscopic dynamics with allowing the next-nearest-neighbor hops results in dramatic changes. Instead 
of the phase separation into high- and low-density regions in the stationary state the system exhibits a 
continuous transition belonging to the Ising universality class for any driving. The relevant features of phase 
diagram are reproduced by an improved mean-field analysis. 

PACS numbers: 02.70.Lq, 05.50.+q, 05.70.Ln, 64.60.Cn 
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The concept of universality is well established in equi- 
librium critical phenomena. According to this concept 
only just a few parameters, i.e. spatial dimension of 
the system and the dimensionality of the order param- 
eter determine the critical exponents meanwhile other 
details like the microscopic dynamics are irrelevant. The 
nonequilibrium systems exhibit richer and more complex 
feature. One of the important questions to address is 
whether this concept also applies to nonequilibrium sys- 
tems. There are examples when a slight extension of 
microscopic dynamics results in different morphology of 
the stationary state . The importance of the dynam- 
ics also manifests in other problems, such as chemically 
reactive mixtures || or driven diffusive systems (DDS) 
jlj where the microscopic and the supposed macroscopic 
model yield different morphologies. The former nonequi- 
librium system also exemplifies that the resulting sta- 
tionary state may differ significantly from that of the 
corresponding equilibrium model j^]. 

Very recently Dickman has introduced a simple driven 
lattice gas model whith hard-core interaction between 
the particles which excludes the simultaneous occupa- 
tion of the nearest neighbor sites on a square lattice || . 
In the absence of driving this model is equivalent to a 
thoroughly investigated equilibrium model discussed in 
connection with the theory of melting |7|-[l0||. For this 
particular interaction the only (control) parameter of the 
model is the particle concentration. The driven version 
of this nearest neighbor exclusion (NNE) model can be 
related to some traffic and granular flow models as de- 
tailed by Dickman |6|. Despite its simplicity the driven 
NNE model shows remarkable nonequilibrium behavior 
such as the coexistence of a low- and high-density phases 
and the current may decrease in response to increased 
drive. The order-disorder phase transition also changes 
from second to first order in the presence of drive. In 
this Brief Report we slightly modify the microscopic dy- 
namics to explore the robustness of this behavior. We 
demonstrate that the coexistence of two phases disap- 
pears and the phase transition remains continuous for all 
drives by allowing the next-nearest-neighbor hops too. A 



more accurate version of the dynamical mean-field theory 
predicts a phase diagram in good qualitative agreement 
with the results of Monte Carlo (MC) simulations. 

We consider two-dimensional driven lattice gases on 
a square lattice with L x L — N sites under periodic 
conditions. The occupation variables <Xj = (1) if the 
lattice site i is empty (occupied) and the concentration 
is defined as p = ^ <Ji/N . The only interaction is to for- 
bid the simultaneous nearest-neighbor occupancy. Dur- 
ing the time evolution a randomly chosen particle can 
hop to one of its empty neighboring sites satisfying the 
condition of NNE. In the basic model studied by Dickman 
the particles can hop only to the nearest neighbor sites 
while in this extended model the next-nearest-neighbor 
hops are also permitted with a probability as defined in 
Fig. [l]. Here, the value of P is varied from 0.5 to 1. The 
case P = 0.5 corresponds to the isotropic hopping rate 
characteristic to the equilibrium model and P = 1 rep- 
resents the infinitely strong drive. Further details and 
the review of the equilibrium properties can be found in 
Ref. §. 
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FIG. 1. Schematic illustration of possible hopping rates 
of the nearest-neighbor (NN) and the extended next-nearest 
neighbor (NNN) hopping dynamics. The drive is horizontal 
and the directions are chosen with equal probability [1/4 (1/8) 
for the NN (NNN) hops]. 

In the stationary state of this model the particle 
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distribution is disordered if the concentration is low 
enough. For sufficiently high concentration sublattice 
ordering can be observed, that is the particles form a 
chequerboard-like pattern. In other words, if the square 
lattice is divided into two interpenetrating sublattices (A 
and B) then the particles prefer to stay in one of these 
sublattices. The transition between the ordered and dis- 
ordered states is characterized by an order parameter 



sition remains continuous for any drives according to this 
level of approximation. 
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where pa and pb denotes the concentration of particles 
in sublattice A and B. In the disordered phase $ = 0. 
Conversely, $ = 1 (or —1) if all the particles are posi- 
tioned in the sites of sublattice A (B). When decreas- 
ing the concentration the equilibrium system (P = 0.5) 
undergoes an order-disorder transition which is continu- 
ous and belongs to the Ising universality class |ic|] . In 
the driven system, however, this transition becomes first 
order if only nearest neighbor hops are permitted ||. 
The present dynamics conserves the number of particles 
therefore the first order transition is accompaned with 
the coexistence of high- and low-density phases. 

To explore the phase diagram of the extended model, 
we have performed a dynamical mean- field analysis. The 
most relevant details of this laborious technique are given 
in the Appendix of the work by Dickman || at the levels 
of two- and four-site clusters. Using this method one can 
determine the probability of possible configurations on a 
given set of clusters if the particle distribution is homo- 
geneous. Similarly to the NN hop model, the two- and 
four-site approximations are not capable to describe the 
effect of driving as a consequence of the strong constraints 
of NNE. This means that the results are independent of 
P and equivalent to those obtained by using the tradional 
cluster variation method jll 12 for the equilibrium sys- 
tem (P = 1/2). The predictions for the phase transition 
point are p^P = 0.25 at two-point and p^ p = 0.317 at 
four-point level. 

The failure of dynamical mean-field approximation is 
disappointing since previously this method gave quali- 
tatively good phase diagrams for several nonequlibrium 
models |13| . To overcome this shortage we have per- 
formed a six-point approximation used succesfully for a 
similar driven lattice gas |]l4j| . In this case we determine 
the configuration probabilities on 2 x 3 clusters. Tak- 
ing the compatibility conditions and symmetries into ac- 
count the configuration probabilities are described by 17 
parameters whose value are determined by solving nu- 
merically the corresponding set of equations of motion. 
At P = 1/2 this calculation has reproduced the same so- 
lution obtained at the level of four-point approximation. 
For the driven cases, however, the solutions (for $ = 
and $ > 0) are already affected by the value of P. As 
an example, for infinite drive the transition appears at 
Pc P (P = 1) = 0.3026. The resulting phase diagram is 
shown in Fig. 0. We should mention that the phase tran- 




FIG. 2. Phase diagram in the p- P plane. Solid curve is the 
prediction of six-point mean-field approximation. Symbols 
represent Monte Carlo results. Dashed line is guide to eyes. 
The antiferromagnetic (A), and the paramagnetic (P) phases 
are indicated. 

To check these predictions MC simulations are per- 
formed for L = 20, 40, 80, 160, and 320 under peri- 
odic boundary conditions varying the concentration p 
and drive P. To reach the stationary state we have used 
both homogeneous random and ordered initial states. In 
the former initial state the particles are positioned within 
a strip parallel to the drive. 

In agreement with the expectation the modification 
of dynamics does not influence the stationary state if 
P = 0.5 (equilibrium model). The MC data of &(p) 
functions collapse when comparing the results obtained 
for the NN and NNN dynamics at any system sizes. At 
the same time the evolution toward the stationary states 
becomes faster when NNN hops are allowed. Similar ef- 
fect was also reported in another driven diffusive model 
Jl| as well as in domain growth processes . In contrast, 
a significant difference can be observed in the stationary 
states when the systems are driven. While the separation 
of the high- and low-density phases characterizes the NN 
dynamics Q leading to a jammed "herring-bone" struc- 
ture then the stationary state remains homogeneous for 
any drive and concentration when NNN jumps are also 
allowed. The study of order parameter and density pro- 
files as a function of transversal coordinate also support 
that the ordered phase is homogeneous in the case of ex- 
tended dynamics. As a consequence, the current has no 
size-dependence and varies smoothly with the concentra- 
tion and P. The phase diagram obtained by MC simula- 
tions confirm the qualitative prediction of the dynamical 
mean- field theory. Namely, the critical concentration (p c ) 
decreases with P and the transition remains continuous 
for those drives displayed in FigJ|. 

The classification of the critical behavior for a nonequi- 
librium model requires careful analysis. The difficulties 
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are well demonstrated by the investigation of standard 
model fl6|| whose critical behavior is even controversial 
p7[ . Here, we present a finite-size scaling analysis for 
the infinite strong drive (P = 1) in the case of extended 
(NNN) dynamics. For this purpose we have used the 
scaling form Ea| 



*(L,p)=L-f , ' v H Pr L 1 ^) , 



(2) 



where p r = (p — p c )/pc is the reduced concentration. As- 
suming that the noncqulibrium phase transition belongs 
to the class of the equilibrium model (i.e. Ising exponents 
are supposed), nice data collapse is found (see Figure^). 
Here, pf c (P = 1) = 0.350(5). This value is consis- 
tent with an alternating estimate which comes from the 
analysis of the fourth order cumulant of order parameter 
U L = 1- (* 4 ) L /3($ 2 )i U|. Similar critical behavior is 
found for any P < 1 drive. 
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FIG. 3. The finite-size scaling of order parameter of NNE 
model with NNN hopping for P = 1. The slopes of the asymp- 
tote lines are indicated. 

In summary, the microscopic dynamics has been 
slightly extended in a previously introduced driven dif- 
fusive lattice gas where the only interaction is to for- 
bid the simultaneous nearest-neighbor occupancy. De- 
spite the weak modification the nonequlibrium behavior 
changes significantly. The stationary states are found to 
be homogeneous for any drives and densities if the next- 
nearest-neighbor hops are also permitted. In this case the 
system exhibits a sublattice ordering when the particle 
concentration is increased. The value of critical concen- 
tration decreases with the drive and this behavior can 
be explained by the higher level of dynamical mean-field 
approximation. The numerical results support that the 
system is in the same universality class of the equilibrium 
Ising model. This model excmplefics that the noncqulib- 
rium behavior may be significantly influenced by a weak 
modification of microscopic dynamics. 



National Research Fund (OTKA) under Grant No. F- 
30449 and T-33098. A. S. thanks the HAS for financial 
support by a Bolyai scholarship. 



[1] A. D. Rutenberg and C. Yeung, Phys. Rev. E 60, 2710 
(1999). 

[2] A. Szolnoki, Phys. Rev. E, 60, 2425 (1999). 

[3] S. C. Glotzer, D. Stauffer, and N. Jan, Phys. Rev. Lett. 

72, 4109 (1994); S. Toxvaerd, Phys. Rev. E 53, 3710 

(1996). 

[4] F. J. Alexander, C. A. Laberge, J. L. Lebowitz, and R. 

K. P. Zia, J. Stat. Phys. 82, 1133 (1996). 
[5] R.K.P. Zia, L. B. Shaw, and B. Schmittmann, Physica A 

279, 60 (2000). 
[6] R. Dickman, Phys. Rev. E 64, 016124 (2001). 
[7] C. Domb, Nuova Cim. 9, Suppl. 1, 9 (1958). 
[8] H. N. V. Temperley, Proc. Roy. Soc. (London) 74, 183 

(1959); ibid. 74, 432 (1959). 
[9] D. M. Burley, Proc. Phys. Soc. (London) 75, 262 (1960); 

ibid. 77, 451 (1961). 
[10] D. S. Gaunt and M. E. Fisher, J. Chem. Phys. 43, 2840 

(1965). 

[11] R. Kikuchi, Phys. Rev. 81, 988 (1951). 

[12] G. W. Woodbury, J. Chem. Phys. 47, 270 (1967). 

[13] R. Dickman, Phys. Lett. A 122, 463 (1987); A. Szolnoki, 

J. Phys. A 30, 7791 (1997); G. Szabo and I. Borsos, Phys. 

Rev. E 49, 5900 (1994); A. Szolnoki, Phys. Rev. E 62 

7466 (2000). 

[14] G. Szabo and A. Szolnoki, Phys. Rev. E 53, 2196 (1996). 
[15] H. C. Fogedby and O. G. Mouritsen, Phys. Rev. B 37, 
5962 (1988). 

[16] S. Katz, J. L. Lebowitz, and H. Spohn, Phys. Rev. B 28, 
1655 (1983); J. Stat. Phys. 34, 497 (1984). 

[17] H. K. Janssen and B. Schmittmann, Z. Phys. B: Condens. 
Matter 64, 503 (1986); P. Garrido, F. de los Santos, and 
M. A. Mufioz, Phys. Rev. E 57, 752 (1998). 

[18] K. Binder, Z. Phys. B: Condens. Matter 43, 119 (1981). 



The authors thank Ron Dickman for stimulating dis- 
cussions. This research was supported by the Hungarian 



3 



